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Abstract 

In the hadronic decays of Z°, gluon emission leads to the appearance of the 
longitudinal component of the fragmentation function, F^. Measurement of Fl and 
the transverse component, Ft, could thus provide an insight into the gluon frag- 
mentation function. However, hadronization corrections at low x can be significant. 
Here we present a method of accounting for such corrections, using the Jetset 
event generator as illustration. 

1 Introduction 



Studies of fragmentation functions have always been important, since these distri- 
butions can not be predicted theoretically, but only be measured experimentally, 
and consecutively be described by phenomenological models. Hadronic decays of 
7*/Z° provide a particularly convenient set of events for analysis and interpreta- 
tion. Helicity components of the fragmentation function, measured in such events, 
can be used in various QCD studies, e.g., the extraction of the gluon fragmentation 
function, and the evaluation of a s . However, existing theoretical calculation being 
restricted to the perturbative region, hadronization corrections must be taken into 
account. In what follows, methods for applying such corrections using the Jetset 
event generator [0 as an example, will be discussed. 

Consider the angular distribution in the process e + e~ —* 7* /Z° — > qq in its rest 
frame. Assuming that the final quark and antiquark are not charge-tagged, i.e. that 
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the forward-backward asymmetry is not accessed, the cross section can be written 



as 



da 3 2 q\ 3. 2 



d(costf) 8 (l + cos^)a x + - S in^ L . (1) 

Here or (<tl) is the cross section associated with a transverse (longitudinal) gauge 
boson polarization state with respect to the qq axis, and 8 is the polar angle of a 
particle with respect to the incoming lepton axis. To lowest order, only mass effects 
contribute to a non- vanishing <tl, but only for the vector part of the cross section, 
and there only with a coefficient <tl/ot = 2m^/E^ m . Even for the b quarks this 
gives a negligible <tl contribution at the energies around the Z° peak. Therefore, 
o"L effectively starts in 0(a s ) of perturbation theory, associated with the emission 
of gluons. 

Since partons are not directly observable, one may define a hadron-level analogue 
ofEq.®g, 

d 2 a h 3 , 9 „, do"x 3 9 „ dor . . 

-(f + cos 2 e)-l + -sin 2 ^. (2) 



dxd(cos#) 8 dx 4 dx 

Here x would preferably be associated with the energy fraction taken by a hadron, 
xe = 2E/E cm , so that X E = 2 in each event. Experimentally it is more convenient 
to use the momentum fraction x p . The transverse and longitudinal fragmentation 
functions are defined by a normalization to the total cross section cr tot = cjt + ol [f| , 

P ^ 1 dCT T , s 1 dcr L 

F T (x) = — , F L {x) = — . (3) 

(T to t dx cr t ot dx 

The former is dominated by the fragmentation of quark jets, whereas the latter 
receives a major contribution from gluon fragmentation. Therefore an experimental 
determination of Fl(x) is a first step towards an extraction of the gluon fragmen- 
tation function, alternative to what is offered by more direct methods in three-jet 
events ||. Several experimental Fl(x) studies have also been presented ||. 

A complication is that hadrons are not moving in the direction of their imag- 
ined mother parton. Already in lowest order of perturbation theory, for qq two-jet 
events, nonperturbative hadronization gives an effective p± smearing that induces 
a nonvanishing -Fl(^) even where none is expected. Furthermore, the association 
of a hadron to a single mother parton is not in agreement with our current best 
understanding of the hadronization process, where it is rather the colour field be- 
tween a colour-connected pair of partons (a string piece Q, or a cluster Q) that 
mediates the hadron production. Therefore the structure of smearing effects may 
become rather nontrivial. Obviously, the effects are especially important at small 
x, which is also the region where one would hope to have some sensitivity to the 
gluon fragmentation function. The string picture also casts in doubt the concept 
of a gluon fragmentation function defined from inclusive quantities, since the string 
fragmentation of a parton depends on the angles to other colour-connected partons. 

The issue of hadronization corrections to fragmentation functions was addressed 
in [||, ||]. The emphasis was on the o"l that can be extracted from J 1 FL(x)xdx 
rather than on F\ J (x) itself, however. Therefore we here address how hadronization 
affects Fl(x) (and -Ft(x)). One main conclusion is that a simple smearing approach 
is not sufficient to describe hadronization effects. Thus it appears impossible to 
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Figure 1: Ft(x p ) and Fl(x p ) for corrected (smeared) hadron angles (dashed lines) and for 
uncorrected (non-smeared) hadron angles (solid lines). 



define a completely model- independent, hadronization-smearing-corrected Fl(x), 
that could be used to extract a gluon fragmentation function. We further suggest 
a correction procedure, based on a cluster search strategy, that should give a less 
model dependent F^(x), but at the price of introducing the cluster resolution scale 
y as a new parameter in the problem. 



2 The simple smearing 

The string model description of qq events introduces a Gaussian transverse momen- 
tum smearing of primary hadrons, oc exp(— p\/2a 2 ) d 2 pj_, where a 0.36 GeV 



1C]. Many primary hadrons are unstable and decay further; this distorts the original 
Gaussian spectrum and reduces the average p±. Since decay products have smaller 
pii, momentum parallel to the jet axis, some correlation is also introduced between 
p± and p\\ . Therefore no simple parameterization is proposed, but instead a Monte 
Carlo simulation with Jetset 7.4 Q is used to histogram the amount of angular 
smearing for different x p bins. 

It is now assumed that this smearing should be applied both to quark and gluon 
jets, so that a hadron will not move exactly in the direction of its mother parton. 
There are obvious shortcomings to equating different kinds of jets, like that gluon 
jets have a lower energy and do not contain decays of charm and bottom hadrons, but 
those particular issues only introduce moderate corrections. More severe objections 
can be raised to the association of hadrons to individual partons, as we will discuss 
further in the next section, but forget for the moment. 

Had particles not been smeared in p±, but parallel with their parton of origin, 
the shape of the angular distribution 

F[x p , cos 6) = hl + cos 2 6) F T (x p ) + - sin 2 6 F L (x p ) (4) 
o 4 
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in a bin of x p could be used to extract Ft(x p ) and Fi,(x p ) in that bin. The abovemen- 
tioned smearing will now modify this. The two angular shapes, (3/8) (1 + cos 2 9) 
and (3/4) sin 2 #, both normalized to unity, are therefore convoluted with the x p - 
dependent smearing distributions, characterized by a distribution in the smearing 
angle # sm and an isotropic azimuthal distribution <^ sm . That is, a parton at an angle 
9 P will produce a hadron at an angle # h , where 

cos 9 h = cos 9 P cos 6 sm - sin # p sin 6 sm cos <p sm . (5) 

Data can now be fitted both to the "non-smeared" angular distribution form, 
Eq.(|j), and to the convoluted ("smeared") analogue. The resulting distributions for 
Ft and Fl, obtained using the Jetset generated events, are shown in Fig. |l[ The 
effect on Ft and Fl is visible for hadron momenta below 10% of the beam energy. 
The low momentum region is affected the most, giving Fl values reduced up to one 
order of magnitude. 



3 Objections to the simple smearing 

The above smearing procedure is correct to lowest order in a s , i.e. it describes how 
two-jet events can induce a nonvanishing Fl,(x p ). We know, however, that hadroni- 
zation of three-jet events cannot be described in terms of a simple incoherent sum of 
three q, q and g jets. One example is the string/drag effect |ll|, p^| , i.e. that particle 
production is suppressed in the angular region between the q and q and enhanced 



in the other two regions, well confirmed experimentally [13|. High- momentum had- 
rons still essentially follow the separate parton directions, but low-momentum ones 
are significantly affected. These are the ones where the angular smearing effects 
are large to begin with. It is well-known that the string effect leads to more two- 
jetlike events, e.g. in terms of thrust T, than implied by symmetric smearing J14j ]. 
The reason is to be found in the enhanced production of particles between two 
colour-connected partons that are close in angle, leading to them seemingly being 
even closer, e.g. that the opening angle between the reconstructed jets typically is 
smaller than that between the original partons. 

There is a nontrivial topology dependence on string effects, especially when 
multiple gluon emission is considered. The issue is therefore best studied in an 
event generator, by comparing angular distributions on the parton (i =p below) 
and on the hadron (i =h) level. As a simple measure of the jettiness of events 
we use T = 1 — T, defined on the parton level of each event. The integrated xe 
spectrum is decomposed as 

d^totV 1 /-i dV 

XEUXE = 



dT J Jo dTdx£;d(cos6 l ) 
= ^(l + C os 2 9)FUT) + lsm*9Fi(T) , (6) 

where the x^-weighting ensures a common normalization F^(T) + F^(T) = 2 at 
parton and hadron level (values of F^ and F£ most conveniently are obtained by 
weighting each particle with an appropriate angular factor [Q]). Then AF^(T) = 
F^(T)—F^(T) is a simple measure of the hadronization impact on Fl. This quantity 
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Figure 2: Event-by-event hadronization corrections AF^(T) = F^(T) — F£(T) for dif- 
ferent simulations: solid line corresponds to Jetset 7.4 PS, dashed - Jetset 7.4 ME 
(both using the string fragmentation), and dotted - Jetset 7.4 ME with independent 
fragmentation scheme. 



is shown in Fig. ^, for one realistic simulation and two toy ones, for uu events at 
91.2 GeV. In the realistic parton shower is used to generate multiparton 

configurations, followed by string fragmentation. The shower develops down to a 
cut-off scale Qq ~ 1 GeV, so that also events in the first bin, T < 0.01, can contain 
some gluons. The other two histograms are based on 0(a s ) matrix elements, where 
only 2- and 3-parton configurations are generated, with a cut T > 0.01 on the 
latter to avoid the singularities of the 3-parton matrix element. Thus the first 
bin here represents pure 2-parton events. While one simulation is again based on 
string fragmentation, the other assumes isotropic smearing around the jet axes, 



basically the independent fragmentation scheme of Hoyer et al. [15, [L4J]. (The same 



fragmentation parameters, tuned to the shower model, have been used in all three 
cases. A retuning of parameters for the 0(a s ) simulations would have given a 
larger nonperturbative p± width a to cover for the lack of perturbative gluons, and 
so would have implied even larger fragmentation smearing.) 

The isotropic smearing is, as expected, giving a rather constant hadronization 
correction AFl(T). There is some jump up in going from two to three jets that 
are smeared, followed by a slow but steady drop with T, since the longitudinal 
component itself is increasing in importance with T and therefore gives an increas- 
ing hadronization smearing of the longitudinal component onto the transverse one 
rather than only the other way around. By contrast, the string fragmentation pro- 
vides a much steeper drop of AFl with T, kicking in immediately when going from 
two to three partons, and enhanced in the shower simulation relative to the simpler 
0(a s ) one. At large T the overall hadronization correction can even turn negative. 
Averaging over the T spectrum (with mean value (T) ~ 0.05), we conclude that 
the typical hadronization smearing contribution is only about a third of the naively 
expected one, as obtained from two-parton results. (Qualitatively this agrees with 
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and probably explains a similar observation in Q of smaller-than-expected hadro- 
nization corrections when using herwig |l6| .) That is, if hadronization corrections 
are viewed as a power series in a s , the 0{a s ) term is of opposite sign and almost as 
big as the 0(1) one. 

It should be remembered, however, that this is integrated over all x^, and that 
we have no similar way of addressing results in specific x bins, since the parton 
and hadron x spectra are quite different. Thus it is likely that the F^Xp) derived 
in the previous section is an underestimation, just like an Fl(x p ) found without 
any smearing corrections is likely to be an overestimation, but it appears impos- 
sible to find the "correct" F^(xp) without making detailed assumptions about the 
hadronization process. 



4 Clustering 



Given the problems with the above smearing recipe, we introduce a new strategy, 
based on the clustering approach. In a nutshell, we propose to rotate all hadrons to 
the direction of the cluster they belong to, as an approximate way of removing had- 
ronization smearing effects. Only thereafter is F\ J (x) extracted from this modified 
cos 9 distribution. The strategy is explained further in the following. 

In clustering algorithms, nearby hadrons are combined to form clusters/jets, in a 
way that should reflect the underlying partonic state, to some approximation. The 
combination process is controlled by (at least) one separation parameter, call it y m i n , 
such that the final state contains no pair of clusters closer to each other than that. 
Clustering algorithms Ccin be applied also to ct pctrtonic state, and here 2/min 

provides 

a regularization of soft and collinear divergences in the perturbative cross sections. 
It is then meaningful to calculate the distribution of partons at a factorization 
scale fi 2 = y m i n E 2 m , and define scale-dependent fragmentation functions parame- 
terizing the subsequent soft-perturbative and nonperturbative hadronization. The 
latter should obey standard QCD evolution equations, starting from some unknown 
nonperturbative form at a low reference scale. 

Over the years many cluster algorithms have been proposed [0], each with its 



strengths and weaknesses. In this article we adopt the Durham one [18], which is 
a standard for many perturbative calculations. The distance measure between two 
clusters i and j is 

2min(£?,£f)(l-cos%) 
Vii = ^2 > ( 7 ) 



VIS 



so that y/ylj roughly corresponds to the relative transverse momentum, scaled to 
the total visible energy £^is (= E cm for an ideal detector). 

If we begin by considering a simple qq event, it should reconstruct back to two 
clusters, unless y m in has been chosen very small. Since the momentum of a cluster 
is given by the vector sum of its constituent hadrons, it would resum opposite and 
compensating p± kicks imparted to hadrons in the fragmentation process. The 
cluster direction should therefore be a better measure of the qq axis than that 
provided by the individual hadron momenta. It is the angular distribution of this 
axis that relates back to the polarization character of the 7*/Z° — > qq decay, and 
that we want to be reflected in our extracted F^(x) and F^(x). Therefore it would 
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be an improvement to rotate all hadrons in a cluster to sit along the cluster direction. 
That is, the 9 of a hadron is redefined while its x value is unchanged. 

At this level there is no contradiction with the smearing approach studied earlier. 
Then we smeared the simple partonic angular shapes to arrive at realistic hadronic 
ones to compare with data, now we un-smear the hadronic angles to approach 
the simple partonic distributions. There is one advantage, however: the clustering 
approach is not sensitive to the width of the p± distribution, i.e. the a parameter, 
unlike the smearing procedure. Of course, the p± width still affects the typical error 
between the qq and cluster axes. 

When considering multijet production, the y m \ n choice does become relevant, 
with u 2 = y m i n E 2 m acting as a factorization scale, as noted above. For a large i/min 
all activity is clustered into the two quark jets, and neither the gluon structure nor 
would be probed. For y m - m — > each hadron or parton is a cluster unto itself, 
and we are back at the starting point. So obviously some intermediate scale is to 
be preferred. Given that the typical hadronization p± width is ~ 0.4 GeV, with a 
tail to larger values, one would conclude that clustering up to p± ~ 1 GeV would 
be a sensible minimum to eliminate the bulk of the hadronization p± smearing. 
At the Z° peak this translates into y m j n ^ 0.0001. In the upper end, we want to 
stay with a picture of multiple gluon emission as the norm, i.e. retain Fl(x) as an 
inclusive quantity, in order not to overlap with traditional studies of gluon jets in 
exclusive three-jet events. Since the average number of clusters per event is three 
for y m j n 0.0025, we conclude that 0.0001 ^ y m j n ^ 0.0025 is a reasonable range, 
over which to study a scale-dependent Fl(x,/i). 




Figure 3: Event-by- event hadronization corrections Aa^/cr tot = (a^ — o-^)/a t ot for different 
?/min scales ( Jetset 7.4 PS and string fragmentation) . 

In Fig. ^, the distribution of event-by-event x^-weighted and -integrated had- 
ronization corrections Aa^/a to t = ~~ £r L)/°tot ^ s shown for some different £/min 
scales, for events generated with parton showers and string fragmentation (super- 
scripts h and p stand for hadron and parton level, respectively). We note the 
significant width of these distributions, showing that event-by-event fluctuations in 
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the hadronization process are important and can be of either sign. Even if small by 
comparison, the mean (A<7L/ctot) does show a systematic bias, positive for small 
y m i n and negative for large y m - m - That is, at small y m [ n the hadronization smearing 
wins over the string effects, while it is the other way around for large y m j n — but 
remember that this is only true when averaging over many events. Nevertheless, one 
possible criterion for a good choice of y m \ n would be where the two effects cancel, 
which then gives y m \ n ~ 0.0002, i.e. fi 1.3 GeV. While a sensible reference value, 
one should not take this particular value too seriously, since it is for one specific 
model, and for one specific set of model parameters. Somewhat different parameter 
values, like for the parton shower cut-off Qo « 1 GeV, defining the parton level of 
the events studied, could lead to slightly different "preferred" a values. 
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Figure 4: Hadronization corrections AF^(T) = F^(T) — F£(T) for different y m i n scales 
( Jetset 7.4 PS and string fragmentation) . 

It is important to note that we here have been considering the x-integrated 
quantity. This is of relevance if one e.g. would like to extract an a s from an 
o"l measurement, and so not uninteresting. For the purpose of determining the 
differential x distribution, F-^(x), however, one would have reason to fear that any 
bias could have an x dependence that would not be caught. In the string model, a 
string piece connecting two partons is boosted by an increasing velocity vector as 
the relative opening angle between the partons is decreased, and so the string effects 
spread upwards to larger x values. A warning signal is then that AFl does depend 
quite significantly on T, Fig. ||, i.e. clustering does not reduce the T dependence 
noted in Fig. [2], but mainly shifts the overall level. Since T probes the topology 
of events, we also do expect this topology to reflect itself in an x dependence of 
hadronization corrections. As in the previous studies, this dependence is then likely 
to show up mainly in the lower end of the x range. At larger x, hadrons are rather 
well aligned with the jet axes, so, even with x-weighting, the few particles out there 
give a small contribution to the (AFl). In Fig [|, the relative difference between 
the inclusive Fl(x p ) and Ff uster (x p ), obtained by replacing hadron angles with 
cluster angles for different y m in values, is shown. It is clearly seen, indeed, that 
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the hadronization corrections are only important at low x p , unless y m i n is chosen 
too high. For the transverse component of the fragmentation function, Ft(x p ), 
corrections have the same absolute amplitude but the opposite sign. 
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Figure 5: Ratio AF L /F L = (F L 
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On the up side, the string effect has its perturbative equivalent in the colour 
dipole That is, asymmetries also exist in the production of soft gluons around 
the direction of the harder partons of an event. Such soft parton emission, below the 
cut-off scale Qo ~ 1 GeV we have used, would thus largely fill in the same regions as 
the nonperturbative hadron production, and with the same topology dependence. If 



one takes Local Parton- Hadron Duality [19] seriously, this equivalence should come 
very close. Our proposed strategy, to reset the 6 angle of particles to that of the 
cluster they belong to, would be applicable also to such perturbatively calculated 
parton topologies. 



5 Summary 



The coherence phenomenon f20f] kills the concept of gluon fragmentation functions 
that can be defined independently of the environment they are found in. The 
"hump-backed" shape of inclusive x spectra Jl9| , [2l| ] is an excellent illustration: by 
coherence the multiplication of partons /hadrons at small x is much less than if the 
hard partons could radiate/hadronize independently. The immediate consequence 
is that the expected "softer gluon than quark jets" picture is difficult to test. This 
impacts both on studies of gluon jets directly in identified three-jet events and 
indirectly via F^x). In this article we have illustrated some of these issues for the 
latter observable. 

It appears safe to conclude that a straightforward extraction of F\^{x) from 
hadron angular distributions exaggerates the rate of particles at small x that should 
be attributed to gluon jets, since even the hadronization of pure qq events induces a 
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'false' Fi,(x) by p± smearing. We have also here shown that a symmetric smearing 
around jet axes introduces a bias in the other direction, since it misses important 
string/drag effects that tend to make three-jet events more two-jetlike. In summary, 
there is no model-independent extraction of a unique Fl(x), especially not at small 
x values. 

We therefore propose to introduce a scale-dependent quantity Fi J (x,^ 2 ). Parti- 
cles are clustered, e.g. with the Durham algorithm, and thereafter assigned the 9 
angle of the cluster they belong to, while retaining their x value. Thus fj, 2 sets an 'un- 
smearing' scale, below which p± fluctuations are killed. We find that a fx « 1.3 GeV 
gives opposite and compensating p± smearing and string effects in Jetset simula- 
tions. While the exact number certainly is model-dependent, the order is a sensible 
one, given that the average hadronization p± is of the order of 0.4 GeV. If this then 
sets a reasonably lower limit, an upper one is related to the desire to stay away from 
the region of exclusive two- or three-jet events. Over an intermediate range, one 
could imagine several experimental determinations providing the scale dependence. 
We also remind that, so far, our studies have only been intended for LEP1 energies. 
Coverage of a wider energy range, e.g. at LEP2, introduces s as a further scale of 
the process and allows more differential tests. 
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